Turbulence of capillary waves revisited 
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Kinetic regime of capillary wave turbulence is classically regarded in terms of three- wave inter- 
actions with the exponent of power energy spectrum being v — —7/4 (two-dimensional case). We 
show that a number of assumptions necessary for this regime to occur can not be fulfilled. Four- 
wave interactions of capillary waves should be taken into account instead, which leads to exponents 
v — —13/6 and v — —3/2 for one- and two-dimensional wavevectors correspondingly. It follows that 
for general dispersion functions of decay type, three- wave kinetic regime need not prevail and higher 
order resonances may play a major role. 

PACS numbers: 47.35.Pq, 05.45.-a, 47.27.-i 



I. Introduction. The very foundations of the kinetic 
(or statistical) wave turbulence theory have been laid in 
[19] where kinetic equation for 3- wave interactions of cap- 
illary waves has been written out and its stationary so- 
lution in the form of power energy spectrum was found. 

Kinetic wave turbulence theory is based on a num- 
ber of assumptions, most important of them being as 
follows: I. weak nonlinearity (nonlinearity is small but 
non-zero and denned by a small parameter < e <C 1); 

II. randomness of phases (all waves interact with each 
other stochastically); III. infinite-box limit (L/X — )> 00, 
where L is the size of the system and A is characteristic 
wave length); IV. existence of an inertial interval in the 
wavenumber space (feo,fei), where energy input and dis- 
sipation are balanced; V. locality of interactions in the 
/c-space (only waves with wavelengths of the same order 
do interact; VI. interactions are locally isotropic (no de- 
pendence on direction). Under these assumptions, wave 
kinetic equations have stationary solutions in the form of 
energy power spectra Ek ~ k~ v ', v > 0, [21]. 

In the case when dispersion function depends only on 
one dimensional parameter, say the gravity constant g for 
water surface gravity waves or surface tension constant 
a for capillary waves, one can compute v using dimen- 
sional analysis, without solving the corresponding kinetic 
equation, [3]. E.g. for a direct cascade we have: 

z/ = 2a + d-6 + (5-3a- d)/(N - 1), (1) 

where a is defined by the form of dispersion function 
lj ~ k a , d is the space dimension of the system and TV 
is the minimal number of waves constituting a resonance 
interaction. 

It has been first established in the frame of the model 
of laminated turbulence [7] that energy spectra Ek have 
gaps formed by exact and quasi-resonances (that is, res- 
onances with small enough resonance broadening). This 
yields two distinct layers of turbulence in an arbitrary 
nonlinear wave system - continuous and discrete - and 
their interplay generates three possible wave turbulence 
regimes: kinetic (energy transport over scales, [21]), dis- 
crete (dynamical energy exchange within a small finite 



set of distinct modes, [8]), and mesoscopic (energy trans- 
port over scales is already observable although only a fi- 
nite - fairly big - number of distinct modes should be 
taken into account, [20]). These regimes and their main 
features are shown schematically in Fig. 1. 

As it was mentioned above, for kinetic regime to as- 
sumptions I- VI must hold, some of which are not easily 
verified in laboratory. However the advantage in this 
case is that the knowledge of dispersion function in a 
wave system immediately yields the explicit form of en- 
ergy distribution over the scales. On the other hand, the 
only assumption necessary for discrete regime to occur is 
weak nonlinearity (I). 

In this Letter we take capillary waves as an illustrative 
example for showing that the standard approach " decay- 
type dispersion function => 3-wave kinetic regime" is 
not universal. 

2. Three- wave interactions. Three- wave resonance 
conditions for capillary water waves with dispersion func- 
tion uj = ak 3 / 2 read 
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ki + k 2 = k 3 



(2) 



Case 1. Wavevectors kj € Z d 



have integer coordinates 
Vj — 1,2,3 (e.g. wave interactions in a resonator are 
regarded) and d is arbitrary. In this case (2) has no 
solutions for arbitrary dimension d of the wave vectors, 
[6]. 

Case 2. Wavevectors kj E R 1 have real coordinates 
and d = 1. In this case, 

kl /2 + k 3 2 /2 = (]fei + k 2 ) 3/2 ^^=0orfc 2 = 0, (3) 

and one can see immediately that for all positive kj the 
right hand side of (3) is always greater than its left hand 
side if both kj ^ 0. If k\ = k : k 2 = cfc, with some 
constant 1 < c < 10 (cf. V), absolute resonance width 

A A = \(Ji +uj 2 -oj 3 \ 

= |fc 3 / 2 + c 3 / 2 fc 3 / 2 -[(c + l)fe] 3 / 2 | 

= fe 3 / 2 |l + c 3 / 2 -( c +l) 3 / 2 |>^^ 3 / 2 (4) 
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FIG. 1: Schematic representation of wave turbulent regimes. 



is rapidly growing function of k when k —> oo (cf. III). 
In particular, if fci = fc 2 = fc, A a ~ 0.82fc 3/2 . 

Case 3. Wavevectors kj G R 2 have real coordinates, 
d = 2, and all three wavevectors are collinear. This case 
can obviously be reduced to the previous one by an ap- 
propriate rotation of coordinate axes. 



Case 



Wavevectors k 7 G R 2 are real valued and 



non-collinear. One might argue that if in this case a 
great amount of almost collinear wavevectors form ap- 
proximate triads with small resonance width, we still can 
expect manifestation of 3-wave kinetic regime in labo- 
ratory experiments in the form of power energy spectra 
Ek,3 ~ fc _T//4 . This case has been studied numerically 
and the results are as follows. 

2.1. Resonance width. Absolute resonance width A^ 
explicitly depends on ki and considering if it is "small" 
or "large" the value of k\ should, of course, be taken 
into account. It is intuitively clear that for larger vectors 
larger resonance width is tolerable, and vice versa. Rela- 
tive resonance width Ar, allowing to distinguish between 
various wave turbulent regimes, might be introduced in 
a number of ways, e.g. [8, 20] and others; the problems 
with introducing a general cumulative function An are 
discussed in [9], Ch.6. 

To perform numerical study of solutions of (2), 
for a pair of two-dimensional wave vectors fci = 
(mi, ni), &2 = (m2, n 2 ) we define relative resonance 
width as absolute resonance of proportional pair with 
norm 1, understanding by the norm of a pair of two- 
dimensional vectors that of the corresponding vector in 
R 4 : ||(fci,fc 2 ) 



\/m\ + n 2 + m| + n\ so that 
A fl =|((m?+n?) 3 / 4 + (m2+ ^ )3 /4 _ 
((m 1 +m 2 ) 2 + (n 1 +n 2 ) 2 ) 3 / 4 )| 



(5) 



where m* 



m j /\\(ki,k 2 



and nj 



nj/\\(ki,k2)\\ with 



J = 1,2. 

2.2. Wave numbers, norms and angles. Our first 
series of numerical simulations served to cast a first 
glance at distribution of wavevectors satisfying (2) in the 
fc-space, primarily, if they are distributed evenly over 



the computation domain or concentrated in some re- 
stricted subdomains. Calculations were performed on 
Z d grid fragments —50 < mi, ni, 7712, n 2 < 50 or 
< mi, ni, m 2l n 2 < 100. 

Exact resonances (with Ar = 0) are achieved for pairs 
(fci,0) and (0,^2) only (cf. Case 1), while for all other 
pairs (fci, k 2 ) approximate interactions may take place. 
In the Fig. 2, left panel, all wavevectors ki,k2 with non- 
negative coordinates taking part in approximate inter- 
actions are shown, and their distribution appears to be 
fairly even. However, if one of the wavevectors, say ki, 
has non- negative coordinates, all k2 interacting with such 
(middle panel), are distributed in fc-space quite irregu- 
larly, leaving completely empty the third quadrant and 
the most part of the first quadrant. Irregularity becomes 
even more striking if we consider interacting pairs where 
both ki,k2 have non- negative coordinates (right panel). 
The most part of the domain consists of wavevectors not 
participating in interactions, while interacting vectors are 
contained in narrow triangles along the axes. Moreover, 
a simple check shows that wavevectors from the lower tri- 
angle interact only with vectors from the upper triangle 
and vice versa, are shown in the Fig. 2, left panel. 

To characterize the ratios of norms of interacting vec- 
tors (cf. V) and angles between them (cf. VI), for each 
solution we computed the ratio of the vector norms k\jk 2 
and the corresponding angle (fci k 2 ) (see Fig. 3). It can be 
seen immediately that solution set is highly anisotropic 
- angles between interacting wavevectors all belong to 
the narrow band between 75° and 87°, i.e. interacting 
wavevectors are almost perpendicular. Norms of the in- 
teracting wavevectors can differ by more than 2 orders 
(Fig. 3, left panel) - maximal ratio found in our solu- 
tion set is k\jk 2 = 101.8. For more than 10% of all the 
solutions, k\jk 2 > 10. Restriction of our attention to 
interactions of wavevectors with norms of the same or- 
der makes angle anisotropy even more pronounced (Fig. 3, 
right panel) - all angles now lie between 75° and 81°, i.e. 
the band width becomes twice smaller. Standard aver- 
aging by angles spectra, [20], obviously can not give any 
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FIG. 2: Two-dimensional wavevectors ki,k2 satisfying (2). Left panel: Wavevectors with non-negative coordinates which 
interact with vectors with arbitrary (positive or negative) coordinates. Computation domain —50 < rai,ni < 50. Middle 
panel: All wavevectors interacting with those shown in the previous panel. Same computation domain. Right panel: Both 
wavevectors have non-negative coordinates. Computation domain < mi,ni < 100. Wavevectors from the lower triangle 
interact only with vectors from the upper triangle and vice versa. 
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FIG. 3: Dependence of ratio of interacting vectors' norms 
(longer to shorter) on the angle between vectors. Left panel: 
Complete picture in computation domain < mj,rij < 100. 
Right panel: Zoomed presentation of the initial interval (ra- 
tio < 10) of the left panel. Axes X and Y denote angles (in 
grad) and ratios correspondingly. 



reliable information in this case. 

2.3. Resonance curves. Solution distribution irregu- 
larities demonstrated above have an elegant explanation. 
Indeed, let us notice two simple properties of the reso- 
nance set of wavevectors satisfying (2): 

• if a pair (fci, fe) is a solution, then every (cfci, C&2) 
is also a solution for any cEM 

• if a pair (fei, £2) is a solution, then every ro- 
tated pair (Tfci, Tfc 2 ) is also a solution for any 

TgSC 



Therefore, it is enough to compute all vectors ki pro- 
ducing resonant interactions with some given ki, say 
ki = (0, 1) to obtain a clear view of the whole reso- 
nant interaction set. Indeed, all resonance partners of 
ki = (0, 1) constitute a smooth curve shown in Fig. 4. 
This curve, as a function n(ra), starts with a flat re- 
gion n ~ m 3 / 2 (left panel), then becomes steeper and for 
m — >> 00 has asymptotic n ~ m 1//2 (right panel). Notice 
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FIG. 4: Color online. Resonance curve of vector (0, 1) in k- 
space, for dispersion function uj ~ k 3 ' 2 . Left panel: The 
initial segment of the curve: m< 1 =>■ n ~ m 3//2 . Right 
panel: The overall view of the curve: for m ^> 1 => n ~ 

™V2 




FIG. 5: Color online. Resonance curves in A;-space (schematic 
representation). Left panel: For the vector ki = (0,1) all 
vectors ki lie on the interaction curve shown. Right panel: 
Two interacting vectors lie on each other's resonance curves 
reciprocally. Resonance curve of the rotated vector is shown 
by the dashed line. 



that the two asymptotic regions lie a few magnitudes of 
10 apart and can not be illustratively presented in one 
figure; so we proceed with schematic representation (Fig. 
5). 

The tangent to the curve drawn from (0, 0) gives the 

k 2 with the minimal angle (ki k 2 ) ~ 74.9. We also see 
that the unit vector can interact both with vectors of 
arbitrarily small and arbitrarily large norms k^. Notice 



that both for k<i —> and fc 2 — » oo the angle (ki k 2 ) — >• 
7r/2. Now, any vector k G R can be produced by stretch 
and rotation of our unit vector, and its resonance curve is 
obtained by stretching the curve of the unit vector (with 
the same coefficient) and rotation (by the same angle). 
If two vectors interact resonantly, then each of them lies 
on the resonance curve of another (Fig. 5, right panel). 

We may conclude with confidence that conditions for 
3-wave kinetic regime to occur are decidedly violated. 
This means that energy cascades observable in various 
laboratory experiments must have another phenomenol- 
ogy. 

3. Dynamic energy cascade. A novel model of a 
dynamic energy cascade generation has been proposed in 
[9] and is based on resonance clusters of a special form 
which are due to decay instability. Below we present 
briefly this model and show how to apply it to capillary 
waves. 

The phenomenon of decay instability has been encoun- 
tered in various fields and is known under different names 
- parametric instability in classical mechanics, Suhl in- 
stability of spin waves, Oraevsky-Sagdeev decay insta- 
bility of plasma waves, modulation instability in nonlin- 
ear optics, Benjamin- Feir instability in deep water, etc., 
[12]. All these so very different physical systems can be 
described in the same frame using the notion of 3-wave 
interactions 



LJi + cj 2 = ^3, ki + k 2 = k 3 , 
or 4- wave interactions 

ui + cu 2 = cj 3 + cj 4 , ki + k 2 = k 3 H 



k 4 , 



(6) 



(7) 



where notation ujj = cj(kj) is used. For instance, 
Benjamin- Feir and modulation instabilities are described 
by (7) with 



ui 



UJS + AcJ, co> 2 = co> 3 — Aco>, 



where < Aoj < 1, 



ki + k 2 = 2k 3 , cji + U02 = 2cj 3 , 



(8) 



(9) 



parametric instability is governed by (6) with k 3 = 
which yields ki = — k 2 , cj 3 = 2uji and so on. Notice 
that parametric instability is usually interpreted via the 
Mathieu equation, [11], describing the motion of an elas- 
tic pendulum when initially the spring-like motion dom- 
inates the pendulum-like behavior. However, the mod- 
ulation equations of motion for an elastic pendulum in 
three-dimensional space with coordinates x, y, z coincide 
with the dynamical equations for a resonant triad (6), 
e.g. [9]. 

General criterion of the decay instability formulated 
by Hasselmann in [5] reads as follows. In a 3-wave sys- 
tem cj 3 -mode is unstable, if any infinitesimally small ex- 
citation of it generates two new modes according to (6); 



modes with frequencies uo\ and cj 2 are neutrally stable. 
Further on we use notions of A- and P-modes for unstable 
and stable modes correspondingly (A is for active and P - 
for passive) as in [10]. (In plasma physics they are some- 
times called mother and daughter modes, e.g. [17].) In 
a 4-wave system any one mode is neutrally stable. How- 
ever, a quartet of the form (9) (called Phillips quartets) 
"imitates" a triad with A- mode having frequency 2cj 3 . 
Energy percolation within a cluster of triads or Phillips 
quartets can be explained in terms of AA-, AP- and PP- 
connections within a cluster of arbitrary structure, e.g. 
[9,10]. 

Suppose that in a wave system A-mode is initially ex- 
cited, then two P-modes are generated. If one of these P- 
modes happens to be A-mode in another triad or Phillips 
quartet - two new P-modes will be generated, etc. This 
process will generate a sequence of nonlinear modes with 
frequencies 



^1,1 +^2,1 = OL -^ 3 ,1, 

^1,2 + ^2,2 = Ct • CJ 3 , 2 , UJ2,2 = Ct • Ctfyi, 

^1,3 + ^2,3 = OL • CJ 3 , 3 , Ld2,3 = OL • CJ 3 , 2 , 

^l,n + ^2,n = OL • UJ3,m ^2,n = OL • ^ 3 ,( n _l), 



(10) 



where a = 1 in a 3-wave system and a = 2 in a 4- 
wave system. A simple assumption that at each step the 
same part p of the energy Eq of the excited A-mode, say 
pEo, < p < 1, is transferred to the P-mode provid- 
ing connection with the next triad or Phillips quartet, 
immediately gives the energy spectrum for the dynamic 
cascade of the form 

En = P n E<d, with some constant < p < 1. (11) 

If the part of energy pj varies from step to step, 

E n = (iq =lPj )Eo, 0< Pj <l,Vj = l,2,...,n. (12) 

The number of modes in the dynamic cascade depends on 
the energy Eq of the initially excited mode. The process 
" dies" when the energy of P-modes become too small for 
weakly nonlinear processes to start, i.e. IV. is violated. 
This means in particular that for smaller E$ , cascade dies 
out faster. For bigger Eq, dynamical cascade is formed 
only at some initial stage after which periodic or chaotic 
energy exchange among all modes within a cluster will 
be observed, though phases are always coherent, see e.g. 

[2]- 
Dissipation can also be included in (11) as follows: 



E n = (n? =1 ^y)£ or (n« = w:)£ , 



(13) 



where < fij < 1, Vj = 1,2, ..,n, and fijEj is the part 
of energy which is lost due to dissipation at the j-th step 
of the dynamic cascade. 



If initially A-mode with frequency a • uo^^n is excited 
and all P-modes in (10) have non-zero frequencies, an 
inverse dynamic cascade is generated: 



• cj 3 , n > a • cj 3> ( n _i) > ... > a- cj 3> i, 



(14) 



while at each cascade step an A-mode decays into two 
P-modes. 

Now suppose that at the initial step of excitation of 
the A-mode with frequency Gcfyi, a narrow zero- frequency 
band < Aujq <C 1 - gains some non-zero energy 
as observed in laboratory experiments, e.g. [13, 16]). 
Then second harmonic co>3,i/2 is generated via 4- wave 
interactions of the form cji } i = u^i/2 — Aco>o, ui2,i — 
^3,1 /2 + Aujq. As soon as we have two excited modes, 
a 3- or 4-wave process may begin in which two P-modes 
generate an A-mode. 

In particular, for the case of vertical parametric exci- 
tation of capillary waves, cascading cluster (10) takes the 
form 



^3,1/2 + ^3,1/2 = 073,1, 
CJ 3 ,i/2 + CJ 3 ,i = 3^3,1/2, 

cj 3 ,i/2 + 3cj 3 ,i/2 = 2cj 3> i, 



(15) 



withco>i 5 j = cj3 5 i/2, Vj 
dynamic cascade: 



1, 2, ..., n, and generates a direct 



cj 3> i < 3cj 3) i/2 < 2cj 3) i < ... < n • cj 3> i/2. 



(16) 



Summarizing, initial excitation of one mode with fre- 
quency oJexit generates an inverse dynamic cascade while 
zero- frequency mode with non-zero energy "triggers" a 
direct dynamic cascade. In the later case, 3- and 4-wave 
interactions of capillary waves take place simultaneously. 
Both types of interactions are characterized by coherent 
phases due to the resonance conditions (6) and (7). 

Indeed, the evidence of strong four-wave coupling in 
nonlinear capillary waves has been identified in [16] by 
computing tricoherence as 

r 2 = \(F 1 F 2 F 3 F; +2 _ 3 )\*/(\F 1 F 2 F 3 \*)(\F 1+2 -. 3 \*), (17) 

where Fj is the Fourier component of the surface eleva- 
tion at the frequency ujj. In general, tricoherence r 2 can 
change from (no phase coupling) to 1 (coherent phases); 
in experiments reported in [16] the level of tricoherence 
t 2 > 0.5 has been observed. 

4. Kinetic cascade vs dynamic cascade. For com- 
paring energy spectra Ek and E n , it is convenient to 
rewrite E n as E n = b~ n Eo with b = 1/p, b > 1. Thus we 
have to compare functions 71 • b~ x and 72 • x~ a , where 
the magnitudes of parameters a, 6, 71,72 are defined by 
the specific wavesystem. As for a, b > 1 
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FIG. 6: Color online. In both panels, function x -1 ' 5 is shown 
by bold black line. Function b~ x is shown by dashed lines of 
various colors for b = 1.4; 1.6 and 2.3 



Ek > E n in the long run. However, for some combina- 
tions of parameters and in some finite domains in /c-space, 
the opposite relation can take place, Ek < E n ; the spec- 
tra Ek and E n might be quite close and even coincide 
for some k (see Fig. 6, left panel). Main characteristics 
allowing to distinguish between kinetic and dynamic cas- 
cades which can be easily observed in experimental data 
are summarized in the Table below. 



Property 


Ek 


E n 


coherent phases 


no 


yes 


dependence on the initial conditions 


no 


yes 


local interactions 


yes 


no 


existence of inertial interval 


yes 


not important 



lim {x a /b x ) = 0, 



(18) 



As some of the initial energy goes into P-modes not 
taking part in connections within the cascading clus- 
ter and zero-frequency mode transfers energy into k- 
spectrum via interactions of the form (8), phases become 
stochastic, 3- wave interactions with big enough resonance 
width appear and kinetic regime can be developed (shown 
in the Fig. 6, right panel). 

This scenario has been confirmed in laboratory experi- 
ments with parametrically excited capillary waves [18], 
the container shaken at frequencies from 0.5 to 3500 
Hz. Energy contained in a zero- frequency band and a 
dynamic cascade are observed; they contain total en- 
ergy Etot of the system at lower forcing. Kinetic cas- 
cade occurs first at frequencies about 220 Hz and its en- 
ergy grows (with increase of the forcing frequency) from 
0.01 • E to t to 0.23 • E totl while energy contained in the 
dynamic cascade decreases from 0.82 • E tot to 0.46 • E tot . 

The understanding of differences between dynamic and 
kinetic cascades is of the utmost importance for correct 
interpretation of the experimental observations. Thus, 
in [1] weak turbulence of capillary waves in Helium has 
been studied and the formation of a local maximum of 
the wave-spectrum near a viscous cut-off was observed 
(under periodic driving force) and correctly attributed to 
the discrete regime (interactions are non-local). On the 



other hand, "in the inertial range dependence of the peak 
amplitudes on frequency is described well by a power law 
function I u ~ uo~ m with the index m w 3.7. This is in 
agreement with the weak turbulence theory which gives 
the value m = 21/6" ([1], p.032001-3). As 21/6=3.5, the 
observed and predicted indexes differ by about 6%. It 
would be worth to check phase coherence in this data in 
order to understand whether this discrepancy is due to 
the available accuracy of measurements or while in fact 
a dynamic cascade is observed and not a kinetic one. 
5. Conclusions. 

• In the system of weakly nonlinear capillary waves 
three substantially different turbulent regimes can occur: 
discrete, kinetic and mesoscopic. 

Discrete regime can be observed in the form of frozen 
turbulence (distinct modes keeping their initial energies, 
energy flux is absent), or as a resonance clustering (a 
set of distinct modes dynamically exchanging energies 
among themselves). Irrotational capillary waves with dis- 
persion function uj ~ k 3 / 2 demonstrate fluxless modes, 
"there is virtually no energy absorption associated with 
high- wave numbers damping in this case" ([15], p. 107). 
On the contrary, rotational capillary waves have disper- 
sion function uj(k) = — § + \/a\k\ 3 + — with constant 
non-zero vorticity Q, and their discrete regime is gov- 
erned by resonance clustering, [4]. The necessary and 
sufficient condition for discrete regime to occur is small 
nonlinear ity, I. 

In kinetic regime energy transport is characterized by 
power spectra of energy (1), it might be observed only if 
the conditions I-VI take place. 

In mesoscopic regime energy transport is realized by a 
dynamic cascade according to (10). 

• Direct dynamic cascade is triggered by 4- wave inter- 
actions of the form (8), i.e. by the presence of a zero- 
frequency mode with non-zero energy. The fact that 4- 
wave resonances play substantial role in the time evolu- 
tion of a system of weakly nonlinear capillary waves is 
confirmed by laboratory measurements in parametrically 
driven capillary waves, [16]. 

• In general, if dispersion function o;(k) has decay type, 
this only means that 3- wave resonance conditions 



o;(ki) + cj(k 2 ) = cj(k 3 ), ki + k 2 = k 3 



(19) 



can be outlined in terms of resonance curves similar to 
those shown in Fig. 5. 
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may have some solutions. This does not necessary mean 
that these solutions possess the properties I-VI. In par- 
ticular, if o;(k) ~ & 7 , 7 > 1, then both properties for- 
mulated in Sec. 2. 3 hold and the geometry of resonances 
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